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ABSTRACT 


Under  the  present  contract  vie  have  developed  a  theory 
for  the  description  and  under standing  of  laser  be¬ 
havior  from  a  fundamental  microscopic  point  of  view. 
Our  model  consists  of  N  two-level  systems  interacting 
with  a  radiation  field.  We  have  treated  problems  of 
collective  motion,  coherence,  dissipation  and  inhomo¬ 
geneous  broadening.  This  report  consists  mainly  of 
the  abstracts  of  six  papers  dealing  with  these  sub¬ 
jects  that  have  or  will  appear  in  scientific  journals 


■tthh 

THEORY  OF  LASER  MODELS 


Introduction 

The  work  performed  under  this  contract  consists  in  developing 
the  consequences  of  a  laser  model  of  N  two-level  systems  in  interaction 
with  a  radiation  field.  The  first  paper  in  the  series  "A  Model  of 
Interacting  Radiation  and  Matter"  introduces  the  Hamiltonian  of  our  model 
and  the  density  matrix  approach  to  the  problem.  In  the  first  paper  we 
treated  atoms  in  the  long  wavelength  limit  and  proved  the  self-consistent- 
field  approximation  is  an  exact  solution  in  a  well-defined  limit.  In 
the  second  paper  Mr.  Picard  and  the  author  of  this  report  developed  the 
coherence  properties  of  our  model.  In  the  third  paper  we  consider  the 
equations  of  motion  of  a  gas  laser  in  the  presence  of  dissipation. 

Paper  four  consists  of  an  extension  of  the  single  mode  case  solved  in 
paper  three  to  the  multimode  gas  laser.  In  paper  five  we  consider  pre¬ 
viously  neglected  pump  induced  center  of  mass  correlations.  Paper  si*: 
extends  the  results  of  paper  four  to  the  case  where  the  average  field 
vanishes.  The  last  paper  is  a  quantum  mechanical  solution  of  the  density 
matrix  equation  for  an  inhomogeneous ly  broadened  gas  laser  to  lowest 
order  in  the  dimensionless  coupling  constant.  The  a.uthor  is  currently 
preparing  a  manuscript  which  includes  a  fully  quantum  mechanical  treat¬ 
ment  of  the  density  matrix  for  a  gas  laser  to  fourth  order  in  the  dir  i- 
sionless  coupling  constant. 

Part  I  of  this  final  report  consists  of  seven  abstracts  of 
scientific  articles,  six  of  which  have  been  published  or  are  in  prer  \, 
and  a  seventh  paper  (paper  number  6)  which  has  appeared  so  far  only 


a  scientific  report. 

Part  II  consists  of  three  abstracts  presented  at  American 
Physical  Society  Meetings.  Two  of  the  abstracts  represent  collabora¬ 
tions  -with  graduate  students,  Mr.  R.  Picard  and  Mr.  S.  T.  Scott,  who 
have  worked  for  short  periods  on  the  present  contract. 

Part  I  —  Published  Articles 

1.  "A  Model  of  Interacting  Radiation  and  Matter” 

Charles  R.  Willis,  J.  Math.  Phys.  5,  1241  (1964) 

ABSTRACT 

We  investigate  the  long-time  behavior  of  a  model  consisting  of 
N  two-level  atoms  in  a  lossless  cavity.  The  Hamiltonian  of  our  system 
contains  the  radiation  oscillators  in  addition  to  the  matter  Hamiltonian 

f 

and  the  usual  j  j  *  Adv  interaction  term.  In  order  to  treat  the  system 
perturbatively,  it  would  be  necessary  to  remove  the  tremendous  degen¬ 
eracy  of  the  system.  Since  this-'ls  prohibitively  difficult,  and  since 
we  are  interested  in  the  long-time  behavior  of  the  system,  we  solve  the 
quantum  mechanical  Liouville  equation  directly  for  a  wide  class  of 
physically  important  initial  distribution  functions.  We  show  the  effec¬ 
tive  expansion  parameter  is  7Ny  where  y  is  a  dimensionless  atomic  dipol< 
moment  and  N  is  the  number  of  atoms.  In  the  lowest  order  we  find  the 
self-consistent-field  approximation.  In  the  next  order,  particle-field 
correlations  appear.  We  explicitly  solve  the  equations  of  motion  for 
particle-field  correlations  in  terms  of  the  average  quantities  that  app 
in  the  self-consistent- field  approximation.  We  show  the  self-consisten 
field  approximation  consists  of  five  first-order  differential  equations. 


Next  we  show  the  equations  of  motion  for  the  density  matrix  of  the  system 
correct  to  order  (yNy)  are  equivalent  to  eight  first-order  differential 
equations.  The  three  additional  equations  are  needed  to  describe  the 
three  second  moments  of  the  density  matrix  of  the  electromagnetic  field 
that  appear  in  second  order.  Our  lowest-order  microscopic  equations  are 
equivalent  to  semiphenomenological  theories  and  our  higher-order  equations 
contain  only  the  measurable  second-order  moments  of  the  electromagnetic 
field  in  addition  to  the  variables  that  appear  in  semiphenomenological 
theories . 

2.  "Coherence  in  a  Model  of  Interacting  Radiation  and  Matter" 

Richard  H.  Picard  and  Charles  R.  Willis 
Phys.  Rev.  139,  A10  (1965) 

-ABSTRACT 

Some  preliminary  results  of  a  theoretical  investigation  of  the 
origin  of  laser  coherence  properties  are  presented.  A  model  consisting 
of  N  two-level  atoms  in  a  perfect  cavity  interacting  through  their  dipole 
moments  with  a  single  mode  of  the  radiation  field  is  adopted,  and  the 
resulting  many-body  problem  is  made  solvable  through  the  introduction  of 
the  self-consistent-field  approximation  (SCFA).  It  is  shown  that  the 
SCFA  is  consistent  with  Glauber's  suggestion  that  an  ideal  laser  might 
be  in  a  "coherent  state",  or  pure  eigenstate  of  the  positive-frequency 
part  of  the  electric-field  operator.  A  conservation  law  is  derived  which 
assures  the  positive  definiteness  of  the  field-density  operator,  and  it  i  ■ 
shown  that  Glauber's  coherent  state  is  that  solution  of  the  SCFA  equation 


which  minimizes  the  corresponding  constant  of  the  motion.  When  the  SCFA 
is  modified  by  the  introduction  of  particle-field  correlations,  the  con¬ 
servation  law  is  destroyed,  and  the  coherent  state  is  no  longer  an  allowed 
solution . 


3.  "A  Model  of  Interacting  Radiation  and  Matter,  II” 

Charles  R.  Willis 

Published  in  J.  Math.  Phys.  1984  (1965) 

ABSTRACT 

We  extend  the  investigation  of  the  long  time  behavior  of  a 
mocel  consisting  of  N  two-level  atoms  interacting  with  radiation  to  include 
dissipation,  pumping,  and  center  of  mass  motion.  We  show  that  when  the 
effective  expansion  parameter  7N7’  is  small,  the  self-consistent-field 
approximation  remains  the  solution  to  lowest  order  in  7N7 .  The  inclusion 
of  the  center  of  mass  motion  introduces  another  dimensionless  parameter  0 
into  the  theory.  We  show  that  when  this  parameter  is  small  the  electro¬ 
magnetic  field  amplitude  varies  slowly  on  the  time  scale  of  the  center  of 
mass  motion.  We  solve  the  equations  of  motion  in  this  slowly  varying 
limit  by  an  extension  of  the  Bogoliuboff-Kryloff  theory  of  quasi -linearity 
to  the  problem  of  time  dependent  integral  kernels.  We  find  the  unique 
stable  stationary  state  and  show  that  in  the  slowly  varying  limit  the 
stationary  state  is  approached  independent  of  initial  conditions.  We 
calculate  the  frequency  shift  to  second  order  in  0.  The  first  order 
frequency  shift  is  the  same  as  that  calculated  by  Lamb.  We  compare  our 
steady  state  solution  with  recent  experiments  with  lasers.  We  include 
the  effect  of  collisions  in  the  steady  state. 


"Model  of  Interacting  Radiation  and  Matter.  HI 


Multimode  Gas  Lasers" 
Charles  R.  Willis 
J.  Math.  Phys.,  in  press 


ABSTRACT 

We  e:  investigation  of  the  long  tine  behavior  of  a 

model  consisting  of  N  two* level  atoms  interacting  with  a  single  electro¬ 
magnetic  cavity  mode  to  include  interaction  with  many  cavity  modes.  We 
show  that,  as  a  consequence  of  the  coupling  between  radiation  modes  pro¬ 
duced  by  spatial  density  variations  of  the  population  inversion,  there  is 
no  strictly  stationary  state  possible  for  multimode  behavior.  However, 
we  obtain  a  stationary  state  by  neglecting  the  rapidly  oscillating  term*. 
The  steady-state  population  inversion  is  then  a  solution  of  an  eigen¬ 
value  problem.  The  eigenvalue  determinant  is  a  function  of  the  number 
of  modes  and  the  coupling  between  modes  in  addition  to  the  usual  depen¬ 
dence  on  frequencies  and  relaxation  times.  We  explicitly  solve  for  the 
unique  eigenvalue  in  special  cases.  The  corresponding  eigenvector  gives 
the  steady-state  mode  intensity  ratios.  The  absolute  values  of  the  steady- 
state  intensities  are  determined  by  the  energy  conservation  equation 
generalised  to  include  pumping  and  dissipation.  We  also  calculate  the 
steady-state  frequency  shifts  for  each  mode.  The  mode  frequency  shifts 
are  practically  independent  of  each  other  and  have  the  eaae  functional 
fora  as  the  single-mode  frequency  shifts. 
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5 .  ’'Models  oi*  Interacting  Radiation  and  Matter  for  Gas  Lasers" 

Charles  R.  Willis 

Proceedings  of  the  Ftysics  of  Quantum  Electronics  Conference, 

to  be  published 
ABSTRACT 

We  investigate  a  model  of  N  two- level  atoms  interacting  with 
radiation  in  the  presence  of  dissipation  and  pumping.  First  we  show  that 
the  self-consi stent-field  approximation  is  a  solution  of  our  equations  to 
lowest  order  in  the  dimensionless  coupling  constant  a.  The  definition  of 
a  is  0^ where  ^  is  the  nuafcer  of  two-level  atoms  per  unit  volume, 
rQ  is  the  classical  electron  radius,  and  >  is  the  wavelength  of  the  radi¬ 
ation.  When  we  consider  the  center  of  mass  motion  ve  get  a  second  dimen¬ 
sionless  constant,  3  ®  a(a)0/<uD )  uhere  u>Q  is  the  atomic  frequency  and  cu^ 
is  the  Doppler  width.  When  0  is  small  t.*e  electromagnetic  field  an^li- 
tude  varies  slowly  on  the  time  scale  of  the  center  of  mass  motion.  We 
find  a  unique  stable  stationary  state  and  show  that  it  is  approached 
rapidly,  independer  of  initial  conditions.  We  find  solutions  of  our 
equations  corresponding  to  sero  »»wrage  field  with  the  same  steady  state 
energy  density  and  cavity  frequency  shift  as  in  the  self-consistent^ 
field  approximation.  When  pump  induced  correlations  are  present,  we  also 
show  using  perturbation  theory  that  the  solution  of  our  equations  corre¬ 
sponding  to  aero  average  field  have  the  came  steady  state  energy  density 
uni  cavity  frequency  shift  as  the  average  field  theories  of  Lamb  n:.i 
Kafc'en  and  raann , 
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’’Models  of  Gas  Lasers" 


Charles  R.  Willis 
To  he  published 
ABSTRACT 

We  extend  the  investigation  of  a  model  of  N  two-level  atoms 
interacting  with  radiation  in  the  presence  of  dissipation  and  pumping 
to  include  physical  situations  where  the  average  fields  vanish  and  where 
purap  induced  correlations  are  important.  The  dimensionless  dynamical 
constant  which  determines  the  nature  of  the  solution  is  p2*  (cu^u^)2where 
is  (Un)  is  the  number  of  two-level  itons  per 

unit  volume,  rQ  is  the  classical  electron  radius,  A  is  the  wavelength  of 
the  radiation  and  tn)Q  is  the  two-level  energy  difference.  The  Doppler 
frequency  characterizes  the  center  of  mass  motion,  in  gas  lasers 
q  «  1  and,  consequently,  the  nonlinear  differential  integral  equations 
which  dewrmine  the  average  value  of  the  electromagnetic  fields  reduce 
to  soluble  nonlinear  differential  equations.  We  show  the  reduction  of 
the  equations  of  motion  apply  to  the  operator  equations  themselves  and 
not  just  to  the  average  values.  Consequently,  we  are  able  to  show  that 
the  steady  state  properties  to  a  given  order  in  9  are  the  same  whether 
th*  average  field  is  zero  or  nonzero*  When  the  pump  induces  correla¬ 
tions  between  the  in  .  .nun!  variables  ami  the  center  of  mass  variables 
of  tho  atcwi,  we  must  use  perturbation  theory  and  our  results  whan  tha 
average  field  if  nonzero  are  equivalent  to  the  theories  of  Lab  sod 
}lm*n  aiid  £&uerMnn.  We  show  that  the  steady  state  properties  with  pump 
induced  correlations  and  average  field  zero  are  equivalent  to  the  results 


•when  the  average  field  does  not  vanish.  Consequently,  steady  state  pheno¬ 
mena  such  as  power  dips,  mode  pulling,  etc.  are  present  when  the  average 
field  vanishes  and  have  the  same  functional  dependence  as  the  nonvanishing 
average  field  results. 

7.  ” Quantum  Theory  of  a  Laser  Model" 

Charles  R.  Willis 
to  be  published 
ABSTRACT 

We  derive  the  kinetic  equations  for  the  coupled  single  particle 

density  matrix  p  and  the  electromagnetic  density  matrix  R  to  lowest  order 

in  the  dimensionless  coupling  constant  p  where  0  =  (o^/cup)  .  The  laser 

frequency  is  (4jt)"^2(^r0A2)^2a>0,  where  ^  is  the  number  of  two- level 

systems  per  unit  volume,  rQ  is  the  classical  electron  radius,  X  is  the 

wavelength  of  the  radiation,  and'Ho^  is  the  two-level  energy  difference. 

The  Doppler  frequency  00^  characterizes  the  center  of  mass  motion.  For 
2 

gas  lasers  0  is  much  less  than  one  and,  consequently,  wc  generalize  and 
use  the  Bogoliubov  derivation  of  kinetic  equations  for  weak  interactions. 
We  find  solutions  when  the  average  field  vanishes  and  which  include  spon¬ 
taneous  emission  correctly.  The  single  particle  density  matrix  ard  the 
radiation  density  matrix  are  coupled  t)  trough  their  second  moments.  When 
we  substitute  the  solution  of  the  sscond  moment  equations  into  the  density 
matrix  equations,  we  find  that  each  density  matrix  satisfies  an  uncoupled 
linear  equation  with  known  time  dependent  coefficients.  We  introduce  and 
uiucuss  dissipation  from  the  density  matrix  point  of  view.  With  the  use 
of  the  density  matrix  formalism  we  indicate  that  the  correct  expansion 
parameter  fer  higher  order  kinetic  equations  is  0>  . 


Fart  II  —  Papers  Presented  at  American  Physical  Society  Meetings 
1.  "Model  of  Interacting  Radiation  and  Matter" 

Charles  R.  Willis 

Bull.  Am.  Pfcya.  Soc .  £,  399  (1964) 

ABSTRACT 


We  have  solved  for  the  lerg-time  behavior  of  a  model  consisting 
of  N  2-level  atoms  in  a  lossless  cavity.  The  Hamiltonian  of  our  system 
contains  the  radiation  oscillators  in  addition  to  the  matter  Hamiltonian 
and  the  usual  J  ^  •  Adv  interaction  term.  In  order  to  treat  the  system 
perturbative  ly,  it  would  be  necessary  to  remove  the  tremendous  degeneracy 
of  the  system.  Since  this  is  prohibitively  difficult  and  since  we  are 
interested  in  the  long-time  behavior  of  the  system,  we  have  solvt-d  the 
quantum-mechanical  Liouville  equation  directly  for  a  wide  class  of  physi¬ 
cally  important  initial-distribution  functions.  We  show  the  effective 
expansion  parameter  is  jUtf,  where  7  is  a  dimensionless  atomic-dipole 
moment  and  W  is  the  number  of  atoms.  In  the  lowest  order,  we  find  the 
self-consistent-field  approximation.  In  the  next  order,  particle-field 
correlations  appear.  We  have  explicitly  solved  the  equations  of  motion 
for  the  particle-field  correlations  in  terms  of  the  average  quantities 
that  appear  in  the  self-consistent-field  aquations. 


2.  "Coherence  in  a  Modal  of  Interacting  Radiation  and  Matter" 

Richard  H.  Picard  and  Char lea  R.  Villi* 

Bull.  Am.  Fbys.  80c.  10,  96  (1965) 

Scam  preliminary  reculte  of  a  theoretical  investigation  of  later 
coherence  properties  era  presented.  A  model  consisting  of  N  2-level 
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systems  in  a  perfect  cavity  interacting  through  their  dipole  moments  with 
a  single  mode  of  the  radiation  field  is  adopted,  and  the  resulting  many- 
hody  problem  is  made  solvable  through  introduction  of  the  self -consistent- 
field  approximation  (SGFA).1  It  is  shown  that  the  3CFA  is  consistent 
with  Glauber's  assumption  that  an  ideal  laser  is  in  a  "coherent  state”, 
or  pure  eigenstate  of  the  positive  frequency  part  of  the  electric -field 
operator.  A  conservation  law  is  derived  that  assures  the  positive  defi¬ 
niteness  of  the  field -density  operator,  and  it  is  shown  that  Glauber's 
state  is  that  solution  of  the  3CFA  equations  that  minimizes  the  corre¬ 
sponding  constant  of  the  motion.  When  the  SCFA  is  modified. by  the  intro¬ 
duction  of  particle-field  correlations,  the  conservation  law  is  destroyed, 
and  the  coherent  state  is  no  longer  an  allowed  solution.  Various  exten¬ 
sions  of  the  results  are  considered  briefly. 

3C.  R.  Willis,  Bull.  Am.  Phys.  Soc.2,  399  (19®0;  J*  Maxh.  Phys.  *>,  12U  (19 64). 
2R.  J.  Glauber,  Phys.  Rev.  1^0,  2^29  (1963) • 

3,  "Dynamical  Coherence  lime  in  a  Gas  Laser” 

Samuel  'r,  Scott  and  Charles  K.  Willie 
Dull.  A?n.  Phys.  See.  10,  1*37  (I'M-). 

We  use  aii  extended  se  If -consistency  argument  to  study  tiynml~ 
cm 1  dephasing  of  the  field  of  a  gas  laser  operating-  in  a  stationary  siwV. 

For  this  purpose,  the  dependence  of  the  atomic- -invars ion  oporatora  or  the 
atomic  velocities  is  orueiaV.  Die  model  includes  dissipation,  piling, 
and  external  radiation  noise.  We  consider  the  atomic  position  operators 


to  be  time-dependent  c  numbers  and  make  the  rotating-wave  approximation 
so  that  the  interaction  Hamiltonian  contains  only  terms  on  the  "energy 
shell".  We  derive  the  equations  of  motion  for  the  atomic-inversion 
operators  and  for  field- correlation  operators  a+(t)a(t-T)  and  average 
over  a  factorized  density  matrix.  We  consider  the  stationary  limit  where 
^a+(t)a(t-T)^>  »  C(t)  and  solve  for  C(t)  self-consistently  in  terms  of 
the  external  pump  and  the  noise  power.  C(t)  determines  the  power  density 
spectrum  and  a  complex  frequency  shift.  The  imaginary  part  of  the  fre¬ 
quency  shift  is  just  the  usual  frequency  shift,  and  the  real  part  is  the 
reciprocal  coherence  time  as  measured  in  a  Miehelson  interferometer.  The 
noise  is  shown  to  contribute  negligibly  to  the  frequency  3hift  and  coher¬ 
ence  time  of  a  gas  laser  above  threshold,  but  causes  a  small  increase  in 
the  pimping  threshold. 
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IS.  AMTNACT 

Wader  the  present  contract  we  have  developed  a  theory  for  the  de¬ 
scription  and  understanding  of  laser  behavior  from  a  fundamental 
microscopic  point  of  view.  Our  model  consists  of  N  two-level  systems 
interacting  with  a  radiation  field.  We  have  treated  problems  of 
collective  motion,  coherence,  dissipation  and  inhomogeneous  broad¬ 
ening.  This  report  consists  mainly  of  the  abstracts  of  six  papers 
dealing  with  these  subjects  that  have  or  will  appear  in  scientific 
Journals. 
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